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1. Introduction 

Sufficiently hot and dense hadronic matter undergoes a transition from a confined hadronic 
phase to a deconfined, chirally symmetric medium, the quark gluon plasma (QGP). At zero chemi- 
cal potential, the transition has been established by lattice QCD simulations, which have proven to 
be a powerful method to analyze the non-perturbative features of the QGP close to the transition re- 
gion. At non-zero chemical potential, however, lattice simulations are limited by the sign problem. 
In fact, direct simulations by standard Monte Carlo methods are not possible. In order to evaluate 
the response of the medium to a small non-zero chemical potential, we perform a Taylor expansion 
of the pressure (grand canonical potential) For the dimensionless combination p/T"^ we define 
the expansion coefficients Cij^ as 

P__P_ _v f ElY f ElV f t!l\'' nn 

74 r%=o ~ ,^/'^H r ; \ t) \t) ' ^ ' 

Here IJ.u,d,s are the chemical potentials of up-, down- and strange-quarks, respectively. The coeffi- 
cients Cij^k have been calculated previously using staggered [Q, ^, ^, ^ and Wilson quarks. As 
they are evaluated at vanishing chemical potential they also provide important information about 
thermodynamic properties of QCD at vanishing baryon number density. 

We present here preliminary results on the smallest non- vanishing coefficients^ which are 
C200 = Cj, coo2 = c|, Clio = Cjf and cioi = c"\.^ Apart from these basic coefficients characterizing 
fluctuations and correlations for different quark flavors, we consider combinations of them that 
yield iso-spin (I) and electric charge (Q) fluctuations, 

ce = i(5c^ + c|-2c^f-cri). (1.2) 

In general, the second order coefficients are related to the quadratic fluctuations, also known as 
susceptibilities of the corresponding charge densities («x)- Fourth order coefficients are related 
to the quartic fluctuations Xa- In terms of the expansion coefficients we have 

^ = (4)-3(4)2 = 24cf with Xe{u,d,s,I,Q,...}. (1.3) 

In Fig. H, we show results for cf and cf for X = m, s calculated with staggered fermions [^]. 
It is apparent that cf varies rapidly and peaks within a small temperature range. This behavior 
may be interpreted as signaling a transition from a confined phase consisting of heavy hadrons to a 

(X) (X) 

deconfined phase consisting of partons. In fact, the ratio xl 1X2 directly related to the squared 
unit charge of quantum number X, for the relevant degrees of freedom [^]. 

The phase diagram of the QGP at /i = is by now believed to be understood quite well. 
E.g., we have plenty of evidence that the transition is a sharp crossover, rather than a true phase 

' It can be shown that the odd-order coefficients vanish through CP-symmetry. 

^In the 2+1 -flavor theory up and down quark masses are degenerate. The coefficients c,jj. and c^/j. are thus equal. 
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Figure 1: Staggered quark results for the Taylor coefficients ^ 
p4fat3 fermions on A^^ = 4 and 6 lattices [0]. 



and c^, obtained with (2+l)-flavor of 



transition. Hence, the general behavior of C2 and C4 as shown in Fig. |l] will not change, even in 
the thermodynamic limit. They may, however be sensitive to the values of the light quark masses. 
At vanishing hght quark mass, the transition at /i = is expected to be second order, in the same 
universality class as the three-dimensional 0(4) -models. In this case the appropriate scaling field 
t (reduced temperature) in the vicinity of the phase transition [||] depends on a combination of T 
and /I 



T-Tr 



Tc 



(1.4) 



This explains the general structure of the temperature dependence of C2„. We thus may expect 
that the behavior of the (2?i)-th jU-derivative of the partition function is similar to that of the n- 
th temperature derivative; e.g. e ~ dlnZ/dT <^ C2 ~ d^lnZ/dji^. A plot of cf versus the 
temperature thus shall resemble the corresponding plot for the energy density, while that of 
shall resemble the specific heat C,,. 

We expect this general pattern of the expansion coefficients also to be reproduced in calcu- 
lations with domain wall fermions. Details of the temperature dependence, in particular at low 
temperature, may, however, be sensitive to chiral properties of the lattice discretization scheme 
used. Our main interest in performing these calculations also with domain wall fermions thus is to 
understand to what extent the better chiral symmetry of this fermion action shows up in thermody- 
namic observables that are directly sensitive to the hadron spectrum. 

1.1 Beyond the Staggered Formulation 

Chiral symmetry plays an important role in determining the nature of the QGP transition. In 
the staggered formulation, the full SU (3)l x SU {3)r chiral symmetry is explicitly broken; at non- 
zero values of the lattice spacing only two U{\) subgroups remain; one in the light and strange 
quark sector, respectively. As a consequence universal properties at the chiral phase transition 
(in the chiral limit) may be recovered only in the continuum limit. Moreover, the loss of full 
chiral symmetry leads to too few light degrees of freedom at low temperatures e.g. one light pion 
instead of three, which may influence bulk thermodynamic properties in this regime. Furthermore, 
questions have been raised about the validity of "the rooting trick" used for staggered fermions [Q], 
which at present are not fully settled |[l^, 1 1 1. 
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It is clearly desirable to go beyond the staggered formulation and work with a fully chiral 
formulation. Domain Wall Fermions (DWF) meet these criteria [[l^]. In this formulation, one 
introduces a fifth dimension {x,y,z,t,s), so that the lattice dimensions are A^^ x x L,. Only 
fermions can propagate along the fifth direction; the gauge fields exist only on the 4-dimensional 
slices 5 = 0, 1, . . .Li — 1. For this reason, the Domain Wall formulation can also be thought of as a 
theory of fermion flavors coupled in a nontrivial way. When one simulates the theory, one finds, 
among other massive modes, two solutions that satisfy [ |T3[ ] 

[Mdn.fiP)]s,' ^J'f ^ (P) = (^-' I sin'P/^j ^J'.v^^ iPl (1-5) 

where M^wf is the Domain- Wall matrix and a sum over s' is implied. Furthermore, these modes are 
eigenstates of 75 and their wavefunctions are found to be localized on the two hyper-planes s = 
and s = Ls — I respectively. These are exactly the chiral modes we want. Using a Pauli-Villars 



subtraction allows to remove contributions of the heavy modes [ ]14[ ] which leaves us with a theory 
that is, for — > 00, chiral-invariant under the full symmetry group. At finite values of a so- 
called residual mass leads to violations of chiral symmetry. In practice one is forced to perform 
calculations at small lattice spacings to suppress these residual mass effects. For thermodynamics 
calculations this means that calculations should be performed for A^t > 8 and/or sufficiently large 
Ls to suppress the residual mass effects. In fact, the calculations we present here only marginally 



satisfy these constraints [ |15[ ] and thus should be considered as a first feasibility study. 

2. Details of the Simulation 

We use here finite temperature 2-i-l-flavor dynamical domain-wall fermion ensembles gener- 



ated by the RBC-CoUaboration [15]. The lattice dimensions are 16^ x 8 x 32. The light and strange 
quark masses used in these calculations are aniq = 0.003 and anis = 0.037, and the domain-wall 
height and anisotropy factor were aMj = 1.8 and a/a^ = 1.00 respectively. The gauge configu- 
rations have been generated at several different values of the coupling, ranging from j8 = 1.95 to 
j8 = 2.14. This covers the transition region from low to high temperatures. On these data sets we 
started to calculate the lowest order expansion coefficients. 

The expansion coefficients can entirely be expressed in terms of traces of the (5-dimensional) 
DWF fermion matrices over space-time, color and spin indices. For the diagonal and off-diagonal 
coefficients we obtain from derivatives of the fermion matrices for light (M„) and strange (M^) 
quarks, 



and 

'-11 



AM\ f I dM 
VT\ \ duxj V dAty, 



cff = 77:;T(tr M-'— tr M-i^ ) ), (2.2) 



respectively. Here we used the shorthand notation M = Mx for X = u, s. We estimated the traces 
using the stochastic random noise method with around 100—150 random vectors for each trace. 
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Products of traces were evaluated in an unbiased manner. Some simulation details and our current 
statistics are summarized in Table 





# conf . 


separation 


#r.v. 


2.14 


40 


10 


100 


2.11 


35 


10 


100 


2.0625 


111 


10 


150 


2.05 


81 


25 


150 


2.0375 


96 


20 


150 


2.025 


71 


20 


150 


2.0125 


125 


10 


150 


1.975 


61 


10 


150 


1.95 


73 


10 


150 



Table 1: Details of the calculation: The columns give from left to right the values of the lattice coupling 
/3, the number of evaluated configurations, the number of trajectories by which these configurations are 
separated and the number of random vectors used for the evaluation of the traces. 



3. Results 

A word about the systematic errors: It can be shown that domain wall fermions satisfy the same 
dispersion relation as naive fermions [|l6|]. This implies that thermodynamic quantities computed 
with the standard DWF action shall have ff{a^) cut-off errors that are of the same magnitude as 
those of a standard staggered or naive fermion discretization scheme. This translates to --^ 10% cut- 
off errors at A^^ = 8 JIt]]. Remarkably, the above continues to be true even after the introduction 



of a chemical potential, with only the magnitude of the conection changing |1 Bp . We note in 
passing that these -dependent cut-off corrections can be computed exactly in the ideal gas limit; 
they are given by Bernoulli polynomials. 

The statistical error in the susceptibilities was determined by the jackknife method. There 
are two independent contributions to the error: The error due to the finite size of the ensemble 
and the error in the trace due to the finite number of random vectors. The latter was found to be 



especially severe for the "disconnected" contributions, which are the products of traces in Eqs. (2.1) 



and (2.2). This shows up in the large error bars for the coefficients Cj, and especially the off- 
diagonal coefficients c "f and c"\ , as can be seen from Fig. ^ Although and do transit from a 
low value to a high one, it is difficult to assign a corresponding transition temperature or equivalent 
a coupling j8 = jSj. 

The disconnected contributions completely or partially cancel each other in the iso-spin and 
electric charge expansion coefficients and C2 ■ This results in much smaller errors for these two 
quantities as can be seen in Fig. |3[ The coefficients and show a smooth transition from a low 
to a high value. There is some indication that they will overshoot the Stefan-Boltzmann (SB) ideal 
gas values at higher temperatures. This would be expected from the analysis of cut-off effects in 



thermodynamics within the DWF formulation [17]. This deserves further analysis. 
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Figure 2: The expansion coefficients Cj, (left) and c"j, c"[ (right) as a function of the coupling j3. 



As we do, at present, have no results on higher order expansion coefficients, e.g. cf, which 
could directly be used to identify a pseudo-critical coupling from the location of a peak in cf , we 



estimate a pseudo-critical coupling, [5 = [5^, from fits of Cj^"''. We use the ansatz [ pj| ] 



JQJ) 



:Atanh{5(/3-i8d)} + C. 



(3.1) 



The resulting curves are also shown in Fig. ^, superimposed on the data points. Table ^ contains 
the best fit values of the free parameters A, B, C and jSj. 

For the transition value of the coupling we obtain = 2.028(6) (jSj = 2.030(7)) from the 
fit to (Cj)- This value is in agreement with the value j3c = 2.031(5) obtained by the RBC- 
collaboration [|^] from an analysis of the chiral condensate. We thus find that the deconfinement 
transition and chiral symmetry restoring transition occur in the same narrow temperature interval. 



From a determination of the Sommer parameter, roa 
able to deduce 



\p=Pj =3.25(18), the group was also 
171(10) (17) MeV [15]. While the latter error arose from the chiral and con- 
tinuum extrapolations, the former was due to the fact that the residual mass varied across the tem- 
perature range and was significantly large for j3 < jSc. 
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A 


B 


C 


A/ 






0.49(5) 


14.1(2.3) 


0.57(3) 


2.028(6) 


1.028 


^2 


0.16(2) 


13.3(2.7) 


0.19(1) 


2.030(7) 


1.418 



Table 2: The best fit values and their errors for the fit of to the ansatz given in Eq. 3.1. 

We have stressed earlier that the current analysis is exploratory and that, in particular, the 
residual mass effects inherent in the DWF ansatz need to be reduced through calculations performed 
for large values of and/or A^^. Clearly, this requires a more thorough analysis in the future. 

Acknowledgments 

We thank Chulwoo Jung and Michael Cheng for helpful discussions. The numerical simu- 
lations were performed on the BlueGene/L computer at the New York Center for Computational 
Science (NYCCS). The work of FK and PH has been supported under Contract No. DE-AC02- 
98CH10886 of the U.S. Department of Energy. 

References 

[1] C. R. Allton et al, Phys. Rev. D 66, 074507 (2002) [arXiv:hep-lat/0204010]. 
[2] C. R. Allton et al, Phys. Rev. D 71, 054508 (2005) [arXiv:hep-lat/0501030]. 
[3] C. Miao and C. Schmidt, Proceedings of this conference. 

[4] R.V. Gavai and S. Gupta, [ArXiv:0806.2233 [hep-lat]]; S. Gupta, Proceedings of this conference. 

[5] C. Bernai-d et al.,Phys. Rev. D 77, 014503 (2008) [arXiv:0710.1330 [hep-lat]]; 
S. Gottlieb, Proceedings of this conference. 

[6] Y. Maezawa et al., PoS LAT2007, 207 (2007) [ai-Xiv:0710.0945 [hep-lat]]. 

[7] S. Ejiri, E Karsch and K. Redlich, Phys. Lett. B 633, 275 (2006) [arXiv:hep-ph/0509051]. 

[8] Y. Hatta and T. Ikeda, Phys. Rev. D 67 (2003), 014028, [ArXiv:hep-ph/0210284]. 

[9] M. Creutz, PoS LAT2006, 208 (2006), [ArXiv:hep-lat/0608020]. 

[10] M. Creutz, PoS LAT2007, 007 (2007) [arXiv:0708.1295 [hep-lat]]. 

[11] A. S. Kronfeld, PoS LAT2007, 016 (2007) [ai-Xiv:0711.0699 [hep-lat]]. 

[12] Y. Shamir, Nucl. Phys. B, 406, 90 (1993), [ArXiv:hep-lat/9303005]. 

[13] G. Fleming, Ph.D. thesis, [ArXiv:hep-lat/0404020]. 

[14] P M. Vranas, Phys. Rev. D 57, 1415 (1998), [ArXiv:hep-lat/9705023]. 

[15] M. Cheng, Proceedings of this conference; D. Renfrew and P Vranas, PoS LAT2007, 235 (2007). 

[16] G. Aarts and J. Foley [UKQCD Collaboration], JHEP 0702, 062 (2007). 

[17] P Hegde, E Karsch, E. Laermann and S. Scheredin, Eur. Phys. J., €55, 423 (2008), 
[ArXiv:0801.4883[hep-lat]]. 

[18] D. Banerjee, R. V. Gavai and S. Sharma, [ArXiv:0803.3925[hep-lat]]. 



7 



